Sensory Attributes and Instrumental Chemical Parameters of Commercial Spanish Cured Ewes’ Milk Cheeses: Insights into Cheese Quality Figures

The external appearance of some of the Protected Designation of Origin (PDO) cured cheeses is similar to other cheese samples made in Spain: 1 kg and 2.5–3 kg formats, cylindrical, and with or without a pleita mark on the surface. In this work, commercial cured ewe’s milk cheese samples with a similar external appearance were analyzed, including five PDO and five non-PDO samples. The parameters analyzed were color, texture, pH, humidity, water activity, and the volatile profile. Additionally, a descriptive and consumer-sensory analysis of the cheese samples was carried out. Statistical analysis of the results showed that luminosity, color coordinates a* and b*, percentage of deformation, humidity, water activity, and acid contents were significantly higher in non-PDO cheese samples. The breaking force, maximum force, and the content of esters were significantly higher in those cheese samples with PDO. In addition, PDO cheese samples showed higher scores for all attributes evaluated by consumers, except for color. These results suggest that PDO cheeses are placed on the market with a higher degree of ripening than non-PDO ones and that consequently they are more positively valued by consumers.


Introduction
The worldwide production of ewe cheese showed an important increase from 645,540 to 683,493 tons between 2010 and 2020 [1].Nowadays, the production of pure ewe cheese in the EU was around 385,444 tons, with Greece being the leading producer (40%) followed by Italy (19%), Spain (18%) and France (16%) [2].Ewe cheese is characterized by a high moisture content (38-43%; p ≤ 0.05) and a lower percentage of crude protein (23-32%; p ≤ 0.01;),in addition to containing a number of total solids around 60% and a value of crude fat ranging from 19 to 29% [3][4][5].
Products with differentiated characteristics due to the origin of their raw materials and/or manufacturing procedures are regulated by European Union Regulation (EC) No. 1151/2012, complemented by Spanish Law 6/201 [6].If these products meet the conditions established by each specific regulation, they may voluntarily adhere to one of the following certifications: protected geographical indication (PGI), protected designation of origin (PDO) and guaranteed traditional specialties (GTS).These three protection schemes favour recognition of the quality and value of the product, mainly by the consumer, and try to combat the unfair competition to which products of a certain reputation are subjected, since the use of protected names for non-protected products not only affects the producer but is considered food fraud for deceiving the consumer [7].Also, these designations constitute a Foods 2024, 13, 127 2 of 16 structural tool in the development and sustainability of rural traditions, as well as a way of promoting Spanish products abroad [8].In this sense, it is reported that PDO cheese samples are sold at higher prices than non-PDO ones [9].
The PDO is a European certification of quality that covers different products whose characteristics and quality are based on the place where they are produced or processed [10].Some of the products covered by the PDO are cheeses, with a total of 26 different cheeses with PDO in Spain.The PDO Regulatory Council allows the production of cheeses in an artisanal way with raw milk and industrially with pasteurized milk.The term "ripened cheese" refers to one that, after the manufacturing process, needs to be kept for a certain time at a suitable temperature and in such conditions that characteristic physical and chemical changes take place.The minimum maturation period is established as 60 days, except for those cheeses made with pasteurized milk and that have a weight ≤ 1.5 kg, for which the minimum maturation is fixed as 30 days with a maximum period of 2 years [11][12][13][14].
The external appearance of PDO cheeses is like other cured ewe cheeses without PDO.However, the final acceptance by consumers depends upon several factors of which flavor is one of the most important [15].This complex attribute depends upon volatile and non-volatile chemical compounds being released from the main fractions of milk during the ripening process: fat, protein, and carbohydrates.Therefore, each product has a characteristic and unique flavor due to its inherent volatile components [16].Thus, this work is focused on the comparative study of instrumental quality parameters (color, texture, pH, water activity, and the volatile-compounds profile) between ewe's milk cheeses with or without PDO.In addition, a descriptive sensory analysis and an affective consumer's study have been carried out to compare the PDO and the non-PDO cheese samples.As a result, the development of fast, sustainable, and reliable techniques for the direct analysis of samples in combination with multivariate analyses is proposed as a useful tool to investigate PDO vs non-PDO cured cheese samples made with ewe's milk, and for the control of cheese samples during food processing to ensure their authenticity.

Cheese Samples
To select the cheese samples for analysis brands with the highest consumption in kg, thus considered the most popular brands, were sampled.That is, the brands sold in the most popular supermarket chains in Spain.The selected cheese samples are all industrially made.Three whole pieces from different batches of those same brands were purchased in single or multiple establishments of those supermarket chains.In this way, differences due to individual pieces of cheese were avoided.So, the studied samples were 30 pieces of cured ewe's milk cheese (3 different batches from 10 different samples) purchased in supermarkets according to the availability offered (Table 1).The following sample inclusion criteria were used: cheese made with pasteurized or raw ewe's milk, minimum maturation time of 30 days, and pieces weighing 3 kg.All the cheese pieces were stored under refrigerated conditions by using vacuum packaging until the day before the analysis.

Sample Preparation
Moisture, volatile compounds, pH, and water activity analyses were run using ground cheese samples.For texture and color analyses, a piece comprising 1/3 of the cheese cylinder was used.Slices 2 cm-thick were obtained after discarding 1.5 cm from the rind.The color was measured from the slices' surface, and immediately afterward the slices were used to obtain cylinders for texture analysis [17,18].For sensory analysis, triangle shape slices were provided to panelists.

Optimization of HS-SPME Procedure
Although HS-SPME analysis of volatile compounds by from ewes' cheeses, such as Manchego, has been previously described, the optimization of this procedure by using a Box-Behnken design (BBD) has not been yet reported [19,20].In this study, three independent factors were considered: extraction temperature (A = 35, 52.5, 70 • C), extraction time (B = 15, 37.5, 60 min), and the amount of sample (C = 0.5, 1.75, 3 g).A total of 17 experiments (3-level design, including 5 center points) were carried out in a randomized order.The sample used for the optimization was a sample of cured ewe's cheese with PDO.The response variable studied was the sum of the total areas of all the volatile compounds extracted, as it was described previously (Table 2).
The design includes the performance of regression analysis of the experimental data so that they fit a second-order polynomial mathematical model described with the following empirical Equation (1): where Y is the predicted response; X represents the system variables; i and j are the design variables; β 0 is a constant; β i is the linear coefficient; β ii is the quadratic coefficient; and β ij is the coefficient of the interactions between variables.

Analysis of Texture, Color, pH, Humidity, and Water Activity
The color was analyzed by using a previously calibrated CM-2600d KONICA MI-NOLTA (Tokyo, Japan) color spectrophotometer, taking 18 random measurements on freshly created surfaces of the cheese, with nine measurements belonging to the outer zone and nine to the interior area.Measurements were taken of the luminosity and the a* and b* color coordinates of the CIELab color space.Cylinders with 1 cm diameter were obtained from the 2 cm-thick slices with a 1 cm diameter punch.In total, 18 cylinders were obtained, nine belonging to the outer zone and nine from the inner zone.All of them were held at 7 • C for 4 h before performance of the uniaxial compression test.The compression of cubes was achieved to 80% of their original height.All tests were carried out at a constant temperature of 7 • C. The texture analyzer used in this study was the TA-TX2 model (Stable Micro System Ltd., Surrey, United Kingdom) with a 245 N load cell and a crosshead speed of 50 mm min −1 .Fracture strength, deformation at fracture, work for fracture, and the maximum force for compression were calculated.
The pH determination was carried out by means of a potentiometric analysis with a Basic 20 Crison Instruments S.A (Barcelona, Spain) model digital pH-meter, using the method proposed by [21].Five grams of cheese was mixed under stirring with 50 mL of distilled water for 5 min to homogenize the sample, and measurements were later taken with the potentiometer electrode.Samples were measured in triplicate.
Moisture was determined by using the method proposed by [21].For this purpose, five grams of cheese were weighed and placed in a JP Selecta Conterm model oven (Barcelona, Spain) for 24 h at 105 • C.After this time had elapsed, the sample was weighed again and the difference in weight was considered the water content of the sample.Samples were measured in triplicate.
The water activity was measured using a previously calibrated AwSprint TH-500 Novasina AG Labware (Lachen, Switzerland) instrument.Cheese samples (grated) were placed in the hygrometer and water activity was measured at 25 • C. Samples were analyzed in duplicate.

Sensory Analysis 2.5.1. Descriptive Sensory Analysis
Ten highly trained panelists (aged 30 to 55 years) with more than 600 h of training in sensory testing from the Department of Agro-Food Technology (UMH) participated in this study (Food Quality and Safety Group).Also, a member of the University of Alicante (Analytical Chemistry, Nutrition and Food Sciences Department) was presented in the two different sessions in which cheeses were evaluated.The selection and training of the panel followed the ISO standard 8586-1 [22].For the present study, the panel worked during 1 orientation session (90 min) during which they trained and discussed the main organoleptic characteristics to be used in the final evaluation.The cheese lexicon used in this session was previously described by other authors [23][24][25][26][27].A unipolar and numerical scale was used to obtain the sensorial information of panelists in which 0 was no intensity and 10 expressed an extremely high intensity.Punctuation of the scale used an increment of 0.5 [28][29][30][31].All samples were presented with 3-digit codes and randomized numbers.
To cleanse the sensations of cheese samples, unsalted crackers and water were provided between samples.

Affective Sensory Analysis
Consumers (n = 100) were recruited among the students/professors of the Food Science and Technology/Agro-Food Engineering Degree at the Miguel Hernández University (UMH, Alicante, Spain) and the Gastronomy Degree at the Alicante University (UA, Alicante, Spain).Their ages ranged from 18 to >45 years old (42%: 18-24 years old; 24%: 25-35 years old; 12%: 36-45 years old; and 22%: >45 years old), being 61% females and 39% males.More than 75% of the consumers were regular consumers of sheep's milk cheese and regular household buyers.An affective consumer's test was carried out following a balanced design of incomplete blocks (split-plot).Authors needed to conduct a study to evaluate 10 products but for this type of evaluation and product type, consumers can give reliable ratings only if 5 products or less are evaluated.Therefore, the described design was selected and, in this way, all of the samples were analyzed by the same number of panelists, but not all the panelists evaluated all the 10 samples.This design was carried out with the support of XLSTAT software [32].The overall satisfaction degree using a 9-point hedonic scale, where 1 meant "dislike extremely" and 9 meant "like extremely", with a "neither like point" denoted 5. Also, questions about the intensity of the attributes using a JAR (just about right) scale were presented.The attributes analyzed were global satisfaction, color, smell, hardness, crumbliness, acidity, saltiness, spiciness, ewe flavor, and aftertaste.

Statistical Analysis
StatGraphics Centurion XV software (version 5.1) was used for the optimization of the HS-SPME procedure (Statistical Graphics Corporation.Rockville, FL, USA).XLSTAT (Microsoft Corporation, Redmond, WA, USA) and StatGraphics Plus 5.0 softwares (Manugistics Inc., Rockville, FL, USA) were used for analyzing the data obtained from the sensory analyses.Statistical tests were performed for the obtained results using the SPSS program with a p-value < 0.05.One-factor ANOVA with a Tukey po-t hoc test were carried out to identify statistically significant differences among the studied samples.

Fitting the Model
Table 2 shows the BBD and the responses based on the sum of total areas of all the volatile compounds extracted from the samples of cured ewe's cheese with PDO.Equation (2) details the polynomial equation of second order, in which the studied response (total area) is expressed as a function of the independent variables (A: extraction temperature; B: extraction time; C: sample weight): Table 3 shows the ANOVA analysis that was carried out to study the reliability of the fitted model.The good precision of the model is shown by the high p-value (0.6067) that indicates no significant ability of the lack-of-fit or the goodness-of-fit of the model to predict the studied response.Also, the R 2 coefficient of the quadratic regression model was 0.7480.

Effects of Independent Variables on Volatiles Extraction
Figure 1 shows that volatiles extraction was only significantly influenced (95% confidence) by the positive effect of the extraction temperature (A).Thus, the extraction yield increased with temperature, and was proportional to the increase in the analyte diffusion coefficient.Bezerra et al. reported similar results when optimization of the volatilecompounds extraction in caprine Coalho cheese was carried out using a response-surface methodology [33].However, a negative A2 interaction effect has been shown.This fact could be related to the partial degradation of some compounds linked to the formation of aromatic compounds due to the Maillard or caramelization reactions.No significant impact on the studied response (p > 0.05) was observed for the rest of the studied parameters.The influence and positive (+) or negative (−) effects of these variables on the response were in the following order: B2 (−) > C (+) > A2 (−) > B (+) > AB (−) > AC (+) > BC (+) > C2 (−).In this sense, individual effects of extraction time (B) and sample weight (C) were positive whereas their quadratic effects (B2 and C2) were negative because of the saturation of the closed vial headspace and, subsequently, the fiber coating [34].Finally, AC and BC interactions were positive, suggesting that a higher sample weight requires higher extraction temperatures and time.However, the negative effect of the AB interaction could promote the degradation of the volatile compounds [35].

Determination of Optimal HS-SPME Conditions
The optimal extracting conditions obtained in the present study were a temperature if 70 °C, a 40 min extraction time, and a 3.0 g sample weight.Under these conditions, the response value predicted by the model was 3.7 × 10 9 .To ensure the reliability of the proposed model, verification experiments under the optimal conditions were performed, in triplicate, obtaining an experimental value of 2.9 × 10 9 ± 0.3 × 10 9 , which did not differ significantly from the predicted value by the model.In conclusion, the developed quadratic model was reliable for the HS-SPME optimization of volatile compounds from hardcured ewe cheeses.

Characterization of the Main Volatile Compounds in Spanish Hard-Cured Ewe Cheeses Using HS-SPME-GC-MS
Firstly, a qualitative study was carried out based on the sum of relative areas to explore differences between samples with and without PDO.More than 55 compounds were identified from the studied samples.Figure 2 shows the average percentage of the classes of volatile compounds found in the cheeses.In general, PDO samples (Figure 2a) were characterized by organic acids (74.2%) followed by ketones (14.7%), aldehydes (7.5%), ethyl (0.9%) and methyl esters (0.7%), and other compounds.In contrast, samples without PDO (Figure 2b) were characterized by a high content of organic acids (91.2%), ketones No significant impact on the studied response (p > 0.05) was observed for the rest of the studied parameters.The influence and positive (+) or negative (−) effects of these variables on the response were in the following order: B2 (−) > C (+) > A2 (−) > B (+) > AB (−) > AC (+) > BC (+) > C2 (−).In this sense, individual effects of extraction time (B) and sample weight (C) were positive whereas their quadratic effects (B2 and C2) were negative because of the saturation of the closed vial headspace and, subsequently, the fiber coating [34].Finally, AC and BC interactions were positive, suggesting that a higher sample weight requires higher extraction temperatures and time.However, the negative effect of the AB interaction could promote the degradation of the volatile compounds [35].

Determination of Optimal HS-SPME Conditions
The optimal extracting conditions obtained in the present study were a temperature if 70 • C, a 40 min extraction time, and a 3.0 g sample weight.Under these conditions, the response value predicted by the model was 3.7 × 10 9 .To ensure the reliability of the proposed model, verification experiments under the optimal conditions were performed, in triplicate, obtaining an experimental value of 2.9 × 10 9 ± 0.3 × 10 9 , which did not differ significantly from the predicted value by the model.In conclusion, the developed quadratic model was reliable for the HS-SPME optimization of volatile compounds from hard-cured ewe cheeses.

Characterization of the Main Volatile Compounds in Spanish Hard-Cured Ewe Cheeses
Using HS-SPME-GC-MS Firstly, a qualitative study was carried out based on the sum of relative areas to explore differences between samples with and without PDO.More than 55 compounds were identified from the studied samples.Figure 2 shows the average percentage of the classes of volatile compounds found in the cheeses.In general, PDO samples (Figure 2a) were characterized by organic acids (74.2%) followed by ketones (14.7%), aldehydes (7.5%), ethyl (0.9%) and methyl esters (0.7%), and other compounds.In contrast, samples without PDO (Figure 2b) were characterized by a high content of organic acids (91.2%), ketones (4.2%), alcohols (1.5%), ethyl esters (1.5%), and other compounds.The most common compounds detected in the samples are detailed in Table 3 and have been previously detected in various cured ewes' cheeses based on previous studies [36][37][38].The ripening process of cheese is very complex and involves microbiological and biochemical changes to the curd that result in the characteristic flavor and texture of each variety.Biochemical changes in cheese samples during ripening may be grouped into primary events (lipolysis; proteolysis; and metabolism of residual lactose, lactate, and citrate) or secondary events (metabolism of fatty acids and of amino acids) [39].Enzymes with lipolytic activity (esterases and lipases) may cause the release of linear-chain acids (hexanoic, nonanoic, octanoic, decanoic and dodecanoic acids) [15].All of them are related to the fatty, cheesy, and waxy flavors of samples.The abundant presence of these short-chain linear fatty acids is thought to be responsible for the characteristic piquant taste of ewes' cheeses [40].Carboxylic acids are also precursors of other aromatic compounds, such as ketones, lactones, aldehydes, and esters [41,42].All of them were extracted and identified in the present work.Consequently, it is expected that PDO samples show a higher ripening stage, due to the high content on secondary aromatic metabolites, in contrast to the non-PDO samples, which are mainly characterized by carboxylic acids.
Using the optimized HS-SPME-GC-MS procedure, 14 volatile compounds were selected for this study based on two criteria.On one hand, the volatile compound had to be present in all the studied samples.On the other hand, identification comprising a minimum similarity of 90% with a compound present in the NIST database was mandatory.Their relative abundance was determined within the selected cheese samples.Identified compounds (Table 4) were grouped under six headings: five carboxylic acids, three esters, three ketones, one aldehyde, one lactone, and one oxime.The ripening process of cheese is very complex and involves microbiological and biochemical changes to the curd that result in the characteristic flavor and texture of each variety.Biochemical changes in cheese samples during ripening may be grouped into primary events (lipolysis; proteolysis; and metabolism of residual lactose, lactate, and citrate) or secondary events (metabolism of fatty acids and of amino acids) [39].Enzymes with lipolytic activity (esterases and lipases) may cause the release of linear-chain acids (hexanoic, nonanoic, octanoic, decanoic and dodecanoic acids) [15].All of them are related to the fatty, cheesy, and waxy flavors of samples.The abundant presence of these shortchain linear fatty acids is thought to be responsible for the characteristic piquant taste of ewes' cheeses [40].Carboxylic acids are also precursors of other aromatic compounds, such as ketones, lactones, aldehydes, and esters [41,42].All of them were extracted and identified in the present work.Consequently, it is expected that PDO samples show a higher ripening stage, due to the high content on secondary aromatic metabolites, in contrast to the non-PDO samples, which are mainly characterized by carboxylic acids.
Using the optimized HS-SPME-GC-MS procedure, 14 volatile compounds were selected for this study based on two criteria.On one hand, the volatile compound had to be present in all the studied samples.On the other hand, identification comprising a minimum similarity of 90% with a compound present in the NIST database was mandatory.Their relative abundance was determined within the selected cheese samples.Identified compounds (Table 4) were grouped under six headings: five carboxylic acids, three esters, three ketones, one aldehyde, one lactone, and one oxime.Table 4. Percentages of peak area of the total volatile composition of cheeses obtained using HS-SPME/GC-MS analysis directly calculated from total ion current (TIC; mean ± standard deviation (SD), n = 3), retention times (min), odor (O), and flavor (F) descriptors between PDO and Non-PDO groups.Three main esters (ethyl hexanoate, ethyl octanoate and ethyl decanoate) were identified in the samples.The presence of esters in cheese samples is related to the esterification of free fatty acids with alcohols both by chemical and enzymatic reactions [44].Ethyl esters are especially known for their important role in the formation of a fruity character in cheese.These compounds can contribute to the aroma of cheese by minimizing the sharpness and the bitterness imparted by fatty acids and amines, respectively [44].The most abundant ester in the studied samples was ethyl decanoate which has been previously described for cured ewes' cheeses [45].The ethyl hexanoate is considered a key odorant in some cheese varieties such as blue cheeses and grana-type cheeses [40] and gives a 'fruity' off-flavor in Cheddar [46], whereas ethyl octanoate seems to be significant in the formation of the aroma of Flor de Guia cheese [44].
Ketones are formed by the enzymatic oxidation of fatty acids to keto-acids and their consequent decarboxylation to methyl ketones [21].In the studied samples, the identified ketones were 2-heptanone, 2-nonanone, and 2-undecanone.2-heptanone was the most abundant ketone detected in this work, in accordance with the literature, which highlighted that this ketone could play an important role in the final fruity, floral, and musty aroma notes of cheeses made from raw ewes' milk such as Zamorano, Manchego and La Serena [19,45].Fruity, floral, and musty notes are associated with 2-nonanone and 2-undecanone, whereas blue cheese notes are attributed to 2-heptanone [44].
Only one aldehyde (phenylacetaldehyde) was identified in all of the studied samples.In effect, aldehydes are transitory compounds that do not accumulate in cheese because they are produced from amino acids either by transamination followed by decarboxylation or by Strecker degradation, being quickly transformed into alcohols or the corresponding acids [44,45].Therefore, a low level of aldehydes indicates a good ripening of the cheese, whereas a high concentration of aldehydes may cause off flavors [42].The aromatic aldehyde phenylacetaldehyde is formed by phenylalanine degradation, via a Strecker reaction [47].This compound has been reported in the literature as the main benzenic compound quantified in Spanish cured ewes' cheeses and has been related to honey-like and floral aromas [44,45].
Only one lactone was detected in all the studied samples, δ-decalactone, which has been found in the highest concentrations in cured ewes' cheeses [45].Lactones have fruity, sweet, creamy, and fermented notes [15].In particular, δ-decalactone contributes to the pleasant coconut odor note [44].δ-Lactones have been reported that may be formed from δ-hydroxyacids following their release from triacylglycerides by lipolysis.So, concentration of lactones usually correlates with lipolysis extent [48].
Finally, there is little information on methoxyphenyl oxime, but the compound has been found in pasteurized cows' milk [49].In a different study, cheeses made with cows' milk with the Streptococcus thermophilus microorganism had a higher content of methoxyphenyl oxime than other cheeses, whereas there was an absence of this compound in directly acidified cheeses, suggesting that it might be produced by bacteria [50].
Based on ANOVA and the Tukey post hoc test (Table 4), it can be confirmed that the volatile profile differs significantly when comparing the studied samples.The variability for all volatiles obtained by Tukey analysis highlighted the necessity of carrying out a multidisciplinary statistical approach in this study.The PCA applied to the 14 volatiles showed two main principal components accounting for 70.37% of the total variation (49.95% F1, 20.43% F2).Projections of sample scores on the plot are shown in Figure 3a.The first function clearly showed differences in the plot space among the samples.In fact, samples with PDO (M1, M3, M4, M9, and M10) were placed closer to each other, showing a homogeneous group in contrast to non-PDO samples (M2, M5, M6, M7, and M8) which were more dispersed in the space.
Taking into account the average values of the identified six chemical families, PDO samples were characterized as having 93.80% composition by acids, 3.72% by oxime compounds, 1.01% by alkyl esters, 0.83% by ketones, 0.46% by lactone and 0.18% by aldehyde compounds, whereas the non-PDO samples showed higher content of acids (94.66%) followed by oxime (2.47%) and ketones (1.57%), with lower contents of esters, lactone, and aldehyde with 0.77%, 0.39% and 0.13%, respectively.All the chemical families showed statistically significant differences when comparing PDO and non-PDO samples (p < 0.05), except for ketones and methoxyphenyl oxime compound.Since carboxylic acids are precursors of other aromatic compounds, such as ketones, lactones, aldehydes, and esters [42], these results seem to indicate that the cheeses with PDO show a more advanced maturation stage, with lower content of carboxylic acids but higher contents of esters, phenylacetaldehyde, and δ-decalactone than the non-PDO samples.This finding agrees with those reported for volatile profiles of Spanish cured cheese with the PDO "Manchego".In this work, samples were studied in two different stages of ripening (4 and 8 months) [51].In conclusion, secondary metabolites such as ketones and esters increased their concentrations during the ripening time.Taking into account the average values of the identified six chemical families, PDO samples were characterized as having 93.80% composition by acids, 3.72% by oxime compounds, 1.01% by alkyl esters, 0.83% by ketones, 0.46% by lactone and 0.18% by aldehyde compounds, whereas the non-PDO samples showed higher content of acids (94.66%) followed by oxime (2.47%) and ketones (1.57%), with lower contents of esters, lactone, and aldehyde with 0.77%, 0.39% and 0.13%, respectively.All the chemical families showed statistically significant differences when comparing PDO and non-PDO samples (p < 0.05), except for ketones and methoxyphenyl oxime compound.Since carboxylic acids are precursors of other aromatic compounds, such as ketones, lactones, aldehydes, and esters [42], these results seem to indicate that the cheeses with PDO show a more advanced maturation stage, with lower content of carboxylic acids but higher contents of esters, phenylacetaldehyde, and δ-decalactone than the non-PDO samples.This finding agrees with those reported for volatile profiles of Spanish cured cheese with the PDO "Manchego".In this work, samples were studied in two different stages of ripening (4 and 8 months) [51].In conclusion, secondary metabolites such as ketones and esters increased their concentrations during the ripening time.

Color, Texture, Humidity, Water Activity and pH Analysis
Texture can be defined as an attribute of cheese originating from a group of physical properties [52].The texture is determined from the breaking force (force needed to break the cylinder, the higher the force the harder the cheese), the deformation at the breaking point (the higher the deformation the more elastic), the breaking work, and the maximum force exerted on the cheeses.A significantly higher mean percentage of deformation of the non-PDO cured cheeses (p < 0.001) was observed, whereas the mean breaking force and the maximum force are significantly higher (p < 0.001) in the cheeses with PDO.In terms of the breaking work, no significant differences were observed between the cured cheeses with or without PDO (Table 5).The obtained results are probably related to the duration of the curing process, since as the maturation advances, the cheese becomes more rigid and less elastic, as has been reported previously by Picon et al. for Manchego cheese [53] and Colin et al. for Cheddar samples in which increased age resulted in increased hardness [54].So, the texture results indicated that the more cured cheeses would be those with the PDO.Texture results were in accordance with those for volatiles, as PDO samples showed lower contents of carboxylic acids but higher contents of esters, phenylacetaldehyde, and δ-decalactone contents than the non-PDO samples.
Maturity of cheese is a complicated process that involves several concurrent and interlinked reactions, including pH changes, lipolysis, and proteolysis [55].Cheese color is also an important physical parameter, which can be considered as an indicator for many attributes like flavor, microbiological quality, and maturity [56].The average color properties of the cheeses, measured using the CIELab color space, show that lightness (L) and the a* and b* coordinates present significantly higher values (p < 0.001) in non-PDO cheeses (Table 5).These differences in color may also be associated with the curing time of the cheese samples, as color is a good indicator of their maturity.The most cured cheeses seem to be those with PDO, as their colors are more yellow-brown as previously described [57].
On the other hand, no significant differences have been detected in the pH values of the analyzed cheese samples whereas the average obtained values of moisture and water activity are significantly higher (p < 0.001) in non-PDO cheese samples as is shown in Table 5.A possible more advanced maturation stage of the cheese with PDO sampled also explains the results as there is a more pronounced loss of water in these samples [58,59].
Therefore, maturation can be considered a process by which the cheese loses water, becoming harder and taking on a stronger flavor.Based on the obtained results, the cheese samples with PDO show, on average, lower values of the color parameters, humidity, and lower water activity, as well as a lower percentage of deformity, but they present a higher breaking force and maximum force.These results seem to indicate that the cheeses with PDO present a more advanced maturation stage than the cured ewe cheeses that do not have the PDO, which is in keeping with volatile compounds results.

Descriptive Sensory Analysis
Thirty attributes showed significant differences among all studied samples.There were statistically significant differences for five visual descriptors (external and internal color, brightness, color homogeneity, and the presence of eyes), six olfactive descriptors (overall intensity PDO ID, lactic, sheep, roasted, fruity (pineapple), and butter), six textural descriptors (touch, firmness, elasticity, crumbliness, graininess, and creaminess), five basic taste descriptors (sourness, saltiness, umami, sweetness, and bitterness), and eight flavor descriptors (intensity, lactic, sheep, roasted, fruity, aftertaste, nutty, and spicy).A partial least squares plot (PLS) was shown to allow a better understanding of the relationships among the ten cured ewes' cheeses, using (i) volatile compounds and (ii) descriptive sensory attributes (Figure 3b).PDO samples (M1, M3, M4, M9, and M10) that can be associated with lactic odor and three flavors (lactic, roasted, and fruity) that can be associated with the three main volatiles situated in the plot near the samples (ethyl hexanoate, phenylacetaldehyde, and δ-decalactone).As was previously described, these compounds are reported in the literature as possible markers of a more advanced maturation stage.Additionally, in Figure 3a, non-PDO samples (M2, M5, M6, M7, and M8) were more dispersed in the space, showing more heterogeneity within this group.

Affective Sensory Analysis
Two PDO samples showed the highest overall similarity values (M3 and M4) whereas the M7 (non-PDO) showed the lowest value (Table 6).M3 also showed the highest value for color with the lowest values shown by M5 and M9.Regarding firmness, again the best satisfaction degrees were obtained for M3 and, in this case, M6 and M10, whereas the lowest ones were those of the samples M2 and M8.The highest crumbliness values were those of the samples M5 and M10, with the lowest values being found for the samples M2, M4, M7, M8, and M9.In general, all samples showed similar aftertaste satisfaction degrees, except for M5 and M6, both non-PDO samples, which showed the lowest values.Sensory profiles of PDO and non-PDO Spanish cured ewes' cheeses reported by the consumer panel are plotted in Figure 4. Significant differences (p < 0.001) between the PDO and non-PDO groups were found for all attributes studied by consumers (Table 6).In general, PDO cheeses showed higher scores than non-PDO ones, except for the color.As a result, it is noticeable that cheeses covered by the PDO obtained a consumer satisfaction degree higher than non-PDO samples.
lowest ones were those of the samples M2 and M8.The highest crumbliness valu those of the samples M5 and M10, with the lowest values being found for the samp M4, M7, M8, and M9.In general, all samples showed similar aftertaste satisfact grees, except for M5 and M6, both non-PDO samples, which showed the lowest va Sensory profiles of PDO and non-PDO Spanish cured ewes' cheeses reported consumer panel are plotted in Figure 4. Significant differences (p < 0.001) between th and non-PDO groups were found for all attributes studied by consumers (Table 6).eral, PDO cheeses showed higher scores than non-PDO ones, except for the colo result, it is noticeable that cheeses covered by the PDO obtained a consumer satis degree higher than non-PDO samples.

Conclusions
HS-SPME-GC-MS, texture, color, humidity, water activity, pH analysis, and sensory analysis have proved to be valuable and reliable parameters for the characterization of different Spanish ewes' cheeses.All of them have been shown to be simple and inexpensive techniques suitable to be used for the control of cheese samples in the food processing industry.The techniques used have shown several advantages such as short experimental time, no use of chemical reagents, and requirement of only a small sample size with minimal sample preparation.
The obtained results suggest that PDO samples had a more controlled production procedure compared to non-PDO samples, as the volatile composition allowed the separation of PDO from non-PDO samples, maybe due to a similar ripening stage of samples.Some indicators were found that could suggest that cheeses with PDO are more mature than other ewe's mature cheeses.For example, lower numbers of carboxylic acids (hexanoic, octanoic, nonanoic, decanoic, and dodecanoic acids) were linked to higher contents of esters (ethyl hexanoate, ethyl octanoate, and ethyl decanoate), δ-decalactone, and phenylacetaldehyde.Also, PDO samples showed higher values for hardness, less water, and more yellowishbrown tones caused by the Maillard reaction, as well as the different textures of these cheese samples that were detected by consumers.In the descriptive analysis, significant differences were observed in 30 attributes.Results obtained from PLS related the PDO samples with lactic, roasted, and fruity flavors, and also with lactic odor, in keeping with the volatile compound results.Attributes analyzed by consumers were better scored for PDO cheeses compared to non-PDO ones.Moreover, the development of the PDO label for Spanish cured ewes' cheeses, showing they have been manufactured and ripened in a defined geographical area, is a valuable tool to help with the preservation and distinction from other industrially made products.

Foods 2024 , 19 Figure 1 .
Figure 1.Pareto chart obtained for total sum of areas of volatiles compounds extracted by HS-SPME.The vertical line indicates the statistical significance at 5% of the effects.

Figure 1 .
Figure 1.Pareto chart obtained for total sum of areas of volatiles compounds extracted by HS-SPME.The vertical line indicates the statistical significance at 5% of the effects.

Foods 2024 , 19 Figure 2 .
Figure 2. Average percentages of the classes of volatile compounds found in PDO (a) and non-PDO (b) cheese samples.

Figure 2 .
Figure 2. Average percentages of the classes of volatile compounds found in PDO (a) and non-PDO (b) cheese samples.

Foods 2024 , 19 Figure 3 .
Figure 3. Multivariate analysis: (a) Plot of the results from Principal Component Analysis performed using the results of the detected volatile compounds; (b) Plot of the PLS obtained by using the results of the relationships between volatiles and odor and flavor descriptors.

Figure 3 .
Figure 3. Multivariate analysis: (a) Plot of the results from Principal Component Analysis performed using the results of the detected volatile compounds; (b) Plot of the PLS obtained by using the results of the relationships between volatiles and odor and flavor descriptors.

Figure 4 .
Figure 4. Sensory profiles of PDO and non-PDO Spanish cured ewes' cheeses obtained by sumer panel.

Figure 4 .
Figure 4. Sensory profiles of PDO and non-PDO Spanish cured ewes' cheeses obtained by the consumer panel.

Table 1 .
Information about the analyzed cured-sheep-cheese samples.

Table 2 .
Combinations of experimental conditions for the BBD and the measured sum of total areas of all the volatile compounds extracted from a sample of cured ewe's cheese with PDO.

Table 3 .
Analysis of variance (ANOVA) for the quadratic model of volatiles extraction.
A, B, and C represent extraction temperature, extraction time, and sample amount, respectively.* Significant p < 0.05.

Table 5 .
Results of the parameters obtained in the texture, color, pH, water activity and humidity analyses.Mean ± Standard Deviation (n = 9).***, significant at p < 0.001.Different superscripts within the same line indicate statistically significantly different values (p < 0.05).

Table 6 .
Consumer satisfaction degree of Spanish cured ewes' cheeses τ , compared between PDO and Non-PDO groups.: not significant at p > 0.05; **, and *** were significant at p < 0.01, and 0.001, respectively.τ Different superscripts for each attribute within the same line indicate statistically significantly different values (p < 0.05). NS